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I. INTRODUCTION 

At first glance topological insulators and quantum spin 
liquids have nothing to do with each other. A topological 
insulator (TI) has properties that can be understood in 
the absence of any electron interactions litSJ while quan- 
tum spin liquids (QSL) are largely understood theoret- 
ically in terms of lattice models of interacting spins 
Moreover, topological insulators explicitly possess charge 
and spin degrees of freedom, while spin liquids are under- 
stood only in terms of spin degrees of freedom. Quantum 
spin liquids come in both gapped and gapless varieties, 
but even for the gapped variety with topological order it 
would appear there is little to say by way of a relation- 
ship between TIs and QSLs. As a supporting argument, 
one might point out that the "type" of topological order 
is different: Topological insulators do not have any non- 
trivial ground state degeneracy, while some of the best 
understood topological spin liquids possess a non-trivial 
ground state degeneracy!^ 

However, to dismiss any relationship between these two 
phases of matter is to ignore some beautiful, and we be- 
lieve insightful, relationships. To help motivate why it is 
useful to look for connections between systems that ap- 
parently have no relationship to one another we need look 
no further than Laughlin's account of his discovery of the 
fractional quantum Hall effect P According to Laughlin, a 
one-dimensional polymer known as polyacetylene (which 
displays charge fractionalization, at least theoretically^) 
was an important motivation for his inspired guess for 
a fractional quantum Hall state wavefunction. 6 ! As it so 
turns out, there are also deep conne ctio ns between poly- 
acetylene and topological insulatorsPS We believe anal- 
ogous relationships between TIs and QSLs can also pro- 
vide useful insights, and could play a key role in deep- 
ening our understanding of strongly interacting variants 
of topological insulators that are not adiabatically con- 
nected to the non- interacting TI limit. 

In this topical review we will flush out some of the 



connections between TIs and QSLs recently established. 
In this paper, we will focus on a number of different lat- 
tice models of fermions. Some of the models are non- 
interacting, some are interacting, some have no intrinsic 
spin-orbit coupling, while some do have intrinsic spin- 
orbit coupling. We study a variety of lattices in two and 
three spatial dimensions and focus on the interplay of 
lattice geometry, interactions, and spin-orbit coupling. 
Our aim is to enlarge our understanding of the phase 
diagram of interacting fermions in two and three dimen- 
sions by solving different model problems and using them 
as points of generalization. For example, which features 
appear often and therefore are rather general, and which 
features appear to be model specific? This approach is 
complementary to directly developing an effective theory. 
Effective theories are extremely powerful because of their 
generality and efficient exploitation of the various sym- 
metries present in a problem. However, it is not always 
straightforward to "guess" an effective theory of a new 
phase of matter so special cases derived from more micro- 
scopic lattice models can serve as a useful guide, partic- 
ularly when strong interactions and/or correlations are 
present. That is the point-of-view we take in this article. 

A summary of our main results is as follows. For 
the case of two-dimensional lattices we show that un- 
der rather general conditi ons th ere is a connection be- 
tween TIs and Kitaev-typ c! 10 * 11 ! spin models. An earlier 
work by Lee, Zhang, and Xiangp] connecting the honey- 
comb Kitaev model and the integer quantum Hall effect 
was an important inspiration for our studies. A relation- 
ship between the integer quantum Hall effect, TIs, and 
Kitaev models exists because the latter can be exactly 
solved by a mapping that reduces the interacting spin 
degrees of freedom to non-interacting Majorana degrees 
of freedom P33 These Majorana fermions hop on the lattice 
and have a band structure similar to "normal" fermions 
(only with a redundancy in the labeling of the states). 
Therefore, one may study the topological properties of 
the energy bands, such as whether they possess a finite 
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Chern number. Energy bands of finite Chern number 
will turn out to establish a connection between TIs and 
Kitaev-type spin liquids (and of course the integer quan- 
tum Hall effect). 

Staying with two-dimensions, we will discuss recent re- 
sults on obtaining a fractional quantum Hall effect in a 
flat-band lattice model with a partially filled band with a 
finite Chern number. This lattice model may prove use- 
ful in the studies of "fractional" topological insulators 
in two-dimensions with time- reversal symmetryP^H^ We 
also describe a class of models in which spin-orbit cou- 
pling can be spontaneously generated at the mean-field 
levelP^ We contrast the relative instability to electron 
interactions of Dirac-like and quadratic band touching 
points in the non-interacting band structure. We find 
that when quadratic band touching points are present, 
topological phases often appear as the leading instabil- 
ity with interactions. Such mean-field studies reduce a 
many-particle problem to a single-particle problem and 
serve at least two key roles in the study of topological in- 
sulators: (1) They elucidate where topological phases can 
be expected to appear in the presence of electron interac- 
tions and (2) They help guide the class of non-interacting 
lattice models that contain interesting physics. We finish 
our discussion of two-dimensional systems with an ex- 
ample of theoretical TI studies guiding the discovery of 
a QSL with novel properties^ and comment on the role 
that disorder can play in stabilizing topological quantum 
spin liquids^ and topological insulators 

In the case of three-dimensions, we focus on sys- 
tems with intrinsic spin-orbit coupling and intermediate 
strength interactions. With an eye towards transition 
metal oxides with Ad a nd 5d elements, we use the slave- 
rotor mean-field theory^EH to investigate the competi- 
tion between interactions, spin-orbit coupling and lattice 
distortions on the pyrochlore lattice. We find that a new 
phase-the "weak topological Mott insulator" -appears^ 
with topologically protected gapless modes that carry 
heat only along a certain class of bulk defects.^ As a 
function of pressure, a complex set of transitions between 
conventional and exotic phases is possible in such sys- 
tems. We discuss the connection of the slave-rotor mean- 
field states to quantum spin liquids describe by fermionic 
excitations and suggest some ways the latter may inform 
the study of interacting TIs in three dimensions.^ We 
close with a discussion of a newly emerging direction 
in transition metal oxides: interaction-driven topologi- 
cal phases at interfacesPBESI 

Our paper is organized as follows. In Sec. |TT] we set 
the convention for our use of the term "topological order" 
in this paper and provide some important general back- 
ground, including TI-QSL connections in one-dimension 
and important results for interacting systems. In Sec. 
Ill we detail our main results for two-dimensions, and in 
Sec. IIVI we describe our chief three-dimensional results. 
Finally, in Sec. [V] we summarize our main results and 
give our view of interesting future directions of study. 



II. PRELIMINARIES 
A. Topological terminology 

Our main goal in this paper is to provide examples of 
how TIs and QSLs are related to each other, and how 
the study of each can help deepen our understanding of 
the other. An important "connection" between these two 
types of phases is via their topological properties. In 
Sec. |III| we will be more precise by what we mean by 
"connection" . We note here that it does not simply mean 
adiabatic connection of the electronic states because TIs 
and QSLs are distinct phases of matter. In order to 
help avoid confusion later, we emphasize that we will use 
"topological" in the broadest sense of the word: there 
is some mathematical invariant characterizing the global 
properties of the system that can not change unless a 
bulk gap in the energy spectrum closes. Those with a 
background in the fractional quantum Hall effect or ex- 
otic quantum spin liquids commonly take "topological 
order" to be synonymous with "non-trivial ground state 
degeneracy" and "non-trivial (i.e. fractional) quantum 
numbers/statistics" of excitations.^ Since a large part of 
this article deals with TIs which do not have a non-trivial 
ground state degeneracy, but nevertheless do have topo- 
logical properties, we choose a more encompassing use of 
the terminology. 

For purposes of this article, we would say that Z 2 time- 
reversal invariant TIs, integer quantum Hall systems, 
fractional quantum Hall systems, and QSLs described at 
low energies by a Z 2 gauge theory all possess "topologi- 
cal order" . We find this use of terminology convenient for 
highlighting the different roles that topological invariants 
and conventional order parameters play in describing a 
state of matter. This way, the quantities characterizing a 
phase are divided into "non-local" (topological) and "lo- 
cal" (conventional order parameter) types. (As an aside, 
we note that such divisions can be subtle and sometimes 
go against conventional wisdom. For example, super- 
conductors should be classified as topologically ordered 
in the sense of WenP as they have a non-trivial ground 
state degeneracy and possess no gauge invariant local or- 
der parameter^*) 



B. General remarks on one-dimensional systems 

In this topical review, we are mainly concerned with 
two and three-dimensional systems. However, a few re- 
marks about one-dimensional systems are in order. One- 
dimensional Fermi systems are special for a number of 
reasons: (i) They are strongly correlated at low-energies 
regardless of the strength of the interactions, and they 
possess interesting properties such as spin-charge sepa- 
ration, a kind of fractionalization of fermion quantum 
numbers! 36 * 37 ^ (ii) Their theoretical study is amenable to 
powerful, non-perturbative techniques bas ed on Abelian 
and non- Abelian bosonization techniques! 36 * 37 ! A well- 
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known gapless phase of one-dimensional fermions conve- 
niently described with these techniques is the Luttinger 
liquid.^ In recent years it has become clear that interest- 
ing, highly universal "spin-incoherent" Luttinger liquid 
regimes also appearPSHUl 

Focusing our attention on gapped one-dimensional sys- 
tems, a few results are suggestive of possible connec- 
tions between TIs and QSLs. While there a re no Zj 
topological insulators in one dimensionpH El th e more 
general Altland and Zirnbauer classifications^*^ 5 -* of non- 
interacting topological states has provided us with im- 
portant results for other topological phases resulting 
from different sets of symmetries. Of particular inter- 
est is the class BDI (time-reversal symmetry with sub- 
lattice symmetry and inte ger sp in), which has a Z clas- 
sification in one dimension! 42 ! 43 ! Fidkowski and Kitaev^l 
have recently shown in a one-dimensional model in class 
BDI that the Z classification breaks down to Z§ in 
the presence of interactions, and have further general- 
ized these results to other symmetry classes.^ In other 
words, phases that are topologically distinct at the non- 
interacting level can be adiabatically connected with in- 
teractions that are slowly turned on and off in a particu- 
lar way. This example proves that interactions can lead 
to important modifications of the non-interacting topo- 
logical classifications! 4 ^ A similar breakdown may occur 
in two and three-dimensional systems, although to our 
knowledge there is not yet a concrete example. 



C. Entanglement as a tool to study topological 
order 

Because topological order is a global, non-local prop- 
erty, it is rather difficult to measure in the general situa- 
tion. (Most physical responses are described by the local 
coupling of an order parameter to an external perturba- 
tion, or generalized "force" P^l) In a few important cases, 
such as the integer and fractional quantum Hall effects, 
and the Z2 topological insulators, the topo logical p roper- 
ties give rise to certain quantized responses! 5 * 8 * 50 * 51 -! These 
examples have topologically protected gapless boundary 
modes that dominate the responses. However, often there 
is no obvious "nice" response that can be computed or 
observed in experiment. The challenge of understanding 
the many-body quantum mechanics in these cases has led 
to important advances in the unde rstanding of quantum 
entanglement in recent yearsP^EH 

In particular, the entanglement entropjEHUl anc j the 
entanglement spectrum^ have emerged as two important 
measures of the quantum entanglement and the topologi- 
cal properties. Typically, a reduced density matrix of the 
ground state is computed by tracing over some degrees 
of freedom (usually real-space coordinates for part of the 
system). The topological entanglement entropy, com- 
puted by tracing over the reduced density matrix times 
the logarithm of the reduced density matrix, contains im- 
portant information about the "quantum dimensions" of 



the excitations in the system.^ 5 -!^! This measure is useful 
for topological systems with excitations that have "ex- 
otic" quantum numbers, such as the fractional quantum 
Hall effect (The non- interacting Z2 topological insu- 
lators do not have "exotic" quantum numbers, so the 
topological entanglement entropy is not a particularly re- 
vealing quantity in this case. However, in a recently iden- 
tified strongly correlated QSH* phase not adiabatically 
connected to the QSH phase, the entanglement entropy 
can serve as a means to distinguish themP^) 

On the other hand, the entanglement spectrum (ES)P^ 
has prov en to be useful in the study of topological 
insulators! 22 * 60 * ^ The ES is found by writing the re- 
duced density matrix as an exponential of a fictitious 
Hamiltonian. The eigenvalues of this fictitious Hamil- 
tonian constitute the entanglement spectrum. It is re- 
markable that the ES, which is determined only from the 
ground-state wave function, contains important informa- 
tion about the excitations of the system.*2Z*In particular, 
for non-interacting models (such as Z2 topological insula- 
tors and superconductors) it has bee n shown that gapless 
boundary modes imply a gapless ES.* 60 * 61 * ^ I Moreover, it 
has recently been proved by Chandran, Hermanns, Reg- 
nault, and Bernevig^ that in fractional quantum Hall 
states the low-energy modes of the entanglement spec- 
trum and edge-state excitation spectrum are in one-to- 
one correspondence. Qi, Katsura, and Ludwig have ob- 
tained similar results.^ 7 -* 

One important result to emerge from the study of the 
ES of topological insulators is that a gapless ES can per- 
sist under some conditions where the physical edge spec- 
trum becomes gapped. For example, applying a magnetic 
field to an inversion symmetric Z2 topological insulator 
will gap the surf ace s tates but leave the entanglement 
spectrum gapless T ne same j s t rue if time-reversal 
symmetry is broken in more subt le wa ysP^ Strong in- 
teractions that gap the edge modea^* 6 ^* will also leave a 
degeneracy in the ESP Thus, the ES "remembers" the 
underlying state is topological. In this case, it turns out 
that the gapless nature of the ES is protected by an in- 
version symmetryP*5Hs31 In the case of three-dimensional 
insulators with inversion symmetry, a gapless entangle- 
ment spectrum implies an "E • B" term in the action 
even if the boundary contains no gapless modest In this 
sense, the entanglement spectrum can be used to identify 
a phase of matter.^ Recently it has been shown that the 
breakdown of the Z classifi cation (to Z$) of class BDI in 
the presence of interaction s] 46 * 47 * can be understood via 
the ESP 

We close this subsection on entanglement by returning 
once more to one dimension and QSLs. For aficianados 
of one-dimensional spin chains, Haldane's classification^! 
of antifcrromagnetically coupled Hcisenberg spins is cele- 
brated: Half-odd integer multiple spin chains are gapless 
(Luttinger liquids) while integer spin chains are gapped 
and possess no long-range order, i.e. they are gapped 
QSLs. Technically, the two cases are distinguished by 
the Berry phase (topological) term in the effective low- 
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energy action. Remarkably, the gapped integer spin 
systems contain topological order which manifests itself 
as a non-local string order parameter with long range 
order anc j protected gapless boundary mode¥^ (simi- 
lar to TIs). In fact, ES studies of integer spin chains 
show phenomenology quite similar to TIs. For example, 
the physical edge state degeneracies lead to the same de- 
generacies in the ES.^ Additionally, perturbations that 
destroy long-range order in the string parameter and 
edge-state degeneracy leave the ES degeneracy in tact, 
so long as inversion symmetry is preserved.^ The ES re- 
sults show that the perturbations to spin chains that de- 
stroy the conventional measures of the topological order 
in the Haldane phase, long-range string order and edge- 
state degeneracies, do not destroy the Haldane phase it- 
self (so long as the bulk gap remains open). (Even in gap- 
less one-dimensional systems the ES may provide useful 
information.^) 

Thus, the entanglement properties have emerged as a 
robust method of identifying a quantum phasePS It is 
remarkable that the behaviors of the ES spectrum for 
TIs and QSLs exhibit such similar phenomenology, both 
with and without perturbations that destroy the edge 
modes. As we will see in the next section, there are more 
direct ways that TIs and QSLs can be related in higher 
spatial dimensions. 



III. TWO-DIMENSIONAL SYSTEMS 

In this section we will describe a class of non- 
interacting and interacting lattice models with topolog- 
ical order in two spatial dimensions. Some of these are 
spin models of QSLs and some are tight-binding models 
of TIs. We will show that under a certain set of con- 
ditions there is a precise topological connection between 
them. 



A. Orientation 

The first prediction of a quantum spin Hall state (a 
two-dimensional TI) in a specific material was made by 
Bernevig, Hughes, and Zhang-- for HgTe quantum wells. 
This prediction was soon verified experimentally in a se- 
ries of beautiful experiments in Wiirzburg! 79 * 80 ! 

In this topical review, we will focus on models ins pired 
by the earlier pioneering work of Kane and MelepHSZl 
which first established the quantum spin Hall state as a 
new topological phase in a simple tight-binding model on 
the honeycomb lattice that preserves time-reversal sym- 
metry. Their Hamiltonian is 

Hk-m = -t X! ( c L c j' ct + H - c -) + ^cdw + Hn 

+ i^S0^2 3tf ■ ^apclaCjp, (1) 



where t is a real first-neighbor hopping, c\ a creates 
a particle of spin a on site i and Ci„ annihilates the 
same particle, Aso is the (real) parameter representing 
the strength of the second-neighbor spin-orbit coupling, 
eij = (djj x df .)/|dy x d?-| is a vector normal to the x — y 
plane describing how the path ((ij)) was tr aversed, that 
is, either "bending" to the right or left! 81 * 82 * Here, "Hcdw 
and "Hr are charge density wave and Rashba spin-orbit 
terms, respectively. The CDW term describes a mean- 
field staggered sublattice potential that would gap the 
Dirac points if Aso = an d put the model in a trivial 
insulator phase. The Rashba spin-orbit term breaks S z 
conservation and is generally present if the physical sys- 
tem has broken inversion symmetry with respect to the 
plane of the honeycomb lattice. 

A key precedent to Kane and Mele's work was a re- 
sult from Haldane^ in which it was shown that a tight- 
binding model of spinless fermions on the honeycomb 
lattice with broken time-reversal symmetry but zero net 
magnetic field could exhibit a quantum Hall effect. (In 
other words, it is the breaking of time-reversal symme- 
try that is key for the quantum Hall effect, not the net 
magnetic flux.) Haldane's model is 

n H = -tJ2(4 c i +H.c.)+^cdw 

(ij) 

+ it 2 ^2 sgn(<% )cjc,, (2) 

((ij)) 

where ti is the (real) parameter describing the string of 
the second neighbor hopping and sgn(eij) is ±1 depend- 
ing on whether the spinless fermion made a right or left 
"turn" along the path ((ij))- 

As Kane and Mele emphasized in their work ISUSSl m 
a quantum spin Hall system that preserves S z quanti- 
zation [H-r. = in ([ij], the state can be represented as 
"two copies" of Haldane's model |2])-one for spin up and 
one for spin down-each of which sees an opposite sign of 
the effective magnetic field and therefore preserves time- 
reversal symmetry overall. The presence of a small but 
finite Rashba coupling does not change the phase; it only 
destroys the S z conservation. 

While the topological classification of the quantum 
Hall effect is based on the Chern number^ (which takes 
the integers Z) and the two-dimensional time-reversal in- 
variant TIs are classified by i^P^ the Haldane-Kanc- 
Mele correspondence shows that the integer quantum 
Hall effect in a lattice model is intimately related to the 
quantum spin Hall state. In particular, if a spin polar- 
ized tight-binding model realizes a quantum Hall effect, 
combining two time-reversed copies of it with spin de- 
grees of freedom is guaranteed to produce a quantum 
spin Hall effect. This simple observation is actually a 
powerful tool for generating desired topological states 
in two-dimensional lattice models. We will use this ap- 
proach in several different ways below, including generat- 
ing TIs from QSLsJS21MI generating QSLs from TIsP (the 
Haldane-Kane-Mele correspondence in reverse) , and gen- 
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erating fractional TIs from the fractional quantum Hall 
effect in lattice models!^ 



B. "Non-interacting" lattice models 

After the work of Kane and Mele on the honeycomb 
lattice^^l there was a flurry of activity exploring simple 
non-interacting tight-binding lattice models on a variety 
of different lattices in two and three spatial dimensions. 
For example, in two dimensions the decorated honeycomb 
lattice,^ the checkerboard lattice,^ the square-octagon 
lattice the kagome lattice,^ the ruby lattice,^ and 
others^ have been shown to support a TI phase in simple 
s-band models. If multiple orbitals on each site are used, 
even the simple square lattice can support a TI.^ 

A natural question to ask is "What properties of a lat- 
tice model would allow it to support a TI phase?" One 
way to view this question in light of the discussion above 
is "What properties of a spinless fermion lattice model 
would allow it to support an integer quantum Hall ef- 
fect?" When cast this way, there are immediate connec- 
tions to so-called "chiral" QSLs, a class of spin liquids 
with broken time-reversal symmetry and a thermal Hall 
conductance. In fact, one of the early models relating 
spin liquids to the (fractional) quantum Hall effect was 
proposed by Laughlin himself!^ 

In recent years, an important class of exactly sol vable 
spin models known as "Kitaev" models has emerged! 10 * 11 ! 
These models can be used to generate concrete examples 
of various types of QSLs (previous studies often relied 
on uncontrolled approximations). In particular, chiral 
spin liquids with a q uantized thermal Hall conductance 
can be foundP^ 90 * 91 * As we now show, such models can 
be adiabatically deformed to an integer quantum Hall 
"Haldane-type" model, which in turn can be used to con- 
struct a two-dimensional TI via the Haldane-Kane-Mele 
correspondence described above. 

As a simple example, consider the Kitaev model on the 
honeycomb lattice,^ 



Kita 



E J<7 K> ( 3 ) 



where the spin-1/2 components <r x ,o~ v , o~ z couple to each 
other along the x, y, z-labeled links, respectively. The 
symmetry of the Hamiltonian ^ is unusual (not Heisen- 
berg, XY, or Ising, for example), but there are recent 
proposals suggestin g it would arise naturally in some 
solid state contexts! 92 * 93 * The main reason for choosing 
the form ([3]) is that it admits an exact solution with a 
spin-liquid ground stated 

The Hamiltonian ^ is easily solved by choosing an 
alternate representation of the spin degrees of freedom: 
One can write the spins as 

(4) 



- r 



where the VjiVjirfjiVj are f° ur "flavors" of Majorana 
fermions on the site j. They have the properties rf = 



iv a V > (v a ) 2 = 1 on each site for a — 0,x,y,z, and they 
satisfy the anticommutation relations {r]f,r]j} = 2S a bSij. 
They contain "half the degrees of freedom of a "normal" 
fermion (one that satisfies the usual anticommutation re- 
lations). Importantly, the transformation Q preserves 
the spin commutation relations [cr°,<T^] = 2it abc a^ pro- 
vided we impose the constraint Dj = ^v'jVjVjVj — 1 
on each site. (The Majorana representation p} actually 
doubles the Hilbert space, so one needs to project to the 
physical Hilbert space.) Under the representation Q the 
Kitaev Hamiltonian is transformed to 



(*.?>» 



'I (5) 



where it; 



is a Z 2 gauge field with eigenvalues 



±1 defined along an a-link. The product of six differ- 
ent Uij around hexagonal plaquettes on the honeycomb 
lattice defines a Zi gauge flux through that plaquette, 
W p = II u ij ■ The usefulness of the representation Q is 
that now 

[Uij,U k l] = [Uij, H K itaov] = [W p , ^Kitaev] = 0, (6) 

so that the u%j are non- dynamical gauge fields and W p 
are non- dynamical gauge fluxes. In other words, the Uij 
in ([5| can be taken to be constants, rather than opera- 
tors where the value of the constants are determined (up 
to gauge equivalent choices) by the fluxes W p around 
the plaquettes in the lattice. For a bipartite lattice (like 
the honeycomb lattice) it is known that the uniform flux 
configuration is lowest in energy,^ and for non-bip artite 
lattices detailed numerical studies indicate the sameP^Hl 

With the Uij as constants, it is clear the interact- 
ing spin Hamiltonian ([3| has now been reduced to a 
non-interacting hopping Hamiltonian |5]). (Recall that 
7 f = (?7 a )^ •) Of course, a non-interacting Hamiltonian 
can always be solved exactly. In the present case, it use- 
ful to Fourier transform to momentum space as one would 
for any tight-binding Hamiltonian and determine the en- 
ergy bands E(k). Because of the self-adjoint property of 
the Majorana fermions, half of the states arc redundant 
and E(k) = —E(-k). One "unusual" property of the 
tight-binding model ([5| is that the hopping amplitudes, 
iJuij, are purely imaginary (assuming J is real). 

The procedure for solving any Kiteav-type model on 
any lattice is exactly the same as above, even if the lat- 
tice is different and/or the spin is larger JHOEDMD Staying 
with the honeycomb lattice model for the moment, it is 
worth emphasizing that Kitaev's solution also determines 
whether the ground state of the spin system sponta- 
neously breaks time-reversal symmetry!^ The statement 
is the following: If the lattice contains plaquettes with 
only an even number of sides (such as the honeycomb lat- 
tice) the ground state will not break time-reversal sym- 
metry. If there are any plaquettes with an odd num- 
ber of links (such as the decorated honeycomb lattice^ 
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or kagome lattice^) the ground state will spontaneously 
break time-reversal symmetry.^ This statement follows 
because time-reversal symmetry sends u.y — > — ity, which 
will change the sign of W p on odd-sided plaquettes, but 
not even-sided ones. 

Applying this to the honeycomb lattice model ^ we 
can immediately infer that since the model does not 
break time-reversal symmetry, the band structure must 
have zero Chern number. On the other hand, the addi- 
tion of the terrrP^l 

H'= J'° V i° Z 3°l+ E J '<rf°H> (7) 

explicitly breaks time-reversal symmetry and leads to a 
second-neighbor hopping term 

U' =^2~iJ'u ij u J kVtVk +^2iJ'u tj u jk T]°7]l, (8) 

{{ik))eA {{ik))ev 

where A and V denote up and down pointing trian- 
gles of (ijk), with i always taken to be the left most 
point. Remarkably, if the Majorana band structure of the 
Hamiltonian ([5| with ^ are taken to represent spinless 
fermions (as opposed to Majorana fermions), the states 
can be adiabatically deformed into the Haldane model ^ 
(which contains purely real first-neighbor hopping and 
imaginary second-neighbor hopping) on the honeycomb 
latticed This implies that the Hamiltonian ^ with |8]) 
is a topological chiral spin liquid with Chern number ±1! 

This result establishes a Kitaev-Haldane-Kane-Mele 
correspondence on the honeycomb lattice that generalizes 
the earlier Haldanc-Kane-Mele correspondence to include 
chiral spin liquids with finite Chern numbers (provided 
the chiral spin liquids can be described by fermions with a 
band structure-there is not such an obvious connected to 
spin liquids described by a bosonic representation of the 
spin degrees of freedom). As we now show, the presence 
of a chiral spin liquid with finite Chern number on the 
honeycomb lattice^ when there is als o an integer quan- 
tum Hall effeclP^ and a TI phase^SH j s a special case of 
a more general result. 

The first exactly solvable model of a chiral spin liquid 
with a finite Chern number that included only two-spin 
interactions was proposed by Yao and Kivelsorp2l on the 
decorated honeycomb lattice shown in Fig.[l] The Hamil- 
tonian is of the Kitaev type and is given by 

>(,u: = E Jafof + E .m;v; + ]T 

(«% {ij)y <«>» 

+ e J'<°* + E JV M + E JV M< ( 9 ) 

(ij)x' 

where the exchange is given by J along the links of the tri- 
angles and J' along the links between triangles, as shown 
in Fig. [I] The Hamiltonian ^ is solved using the trans- 
formation Q and a fermion hopping model with purely 
imaginary hopping parameters results, similar to what 
is found in ([5]). Because of the triangular plaquettes on 




FIG. 1. (color online) The decorated honeycomb (sometimes 
called the "star") lattice. An exactly solvable chiral spin 
liquid model with a finite Chern number and non-Abelian 
anyons has been found on this latticed By the Kitaev- 
Haldane-Kane-Mele correpondence, this implies a topological 
insulator phase in a simple s-band model also exists on this 
latticeP 



the lattice, the ground state spontaneously breaks time- 
reversal symmetry. For < J' / J < y/3, the ground state 
is a gapped chiral spin liquid with Chern number ±1, 
while for J'/ J > the ground state is gapped chi- 
ral spin liquid with zero Chern numberPSl If the Kitaev- 
Haldane-Kane-Mele correspondence described above on 
the honeycomb lattice is to prove general, Yao and Kivel- 
son's discovery of the chiral spin liquid with finite Chern 
number in a Kitaev model should imply there is a TI 
in an s-band model on the decorated honeycomb lattice. 
We have shown this is indeed the case!^ 
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FIG. 2. The decorated honeycomb lattice energy bands 
along high symmetry directions in the Brillouin zone (from 
Ref.[H]). In ( a ) ^so = and in (b) A so = O.lt. There are 
Dirac points at K and K' (not shown) and quadratic band 
crossing points (QBCP) at T in (a), while in (b) Aso 7^ 
opens up a gap at each of these points and destroys the flat 
bands. 
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Consider the s-band Hamiltonian on the decorated 
honeycomb lattice, 

(ij), <r,A (ij) ,<r,A-+A 

(10) 

which describes nearest-neighbor hopping on the trian- 
gles "A" with amplitude t and between triangles "A — » 
A" with amplitude t', and the remaining terms describe 

). The bands for 
2l and the phase 



the same physics as represented in (I 
Hcdw = Wr — are shown in Fig, _ 
diagrams at various filling fractions are shown in Fig. [3] 
Note that at all filling fractions there is a finite region of 
parameter space occupied by a TI (QSH) phase. 



Kitaev-Haldane-Kane-Mele correspondence on the dec- 
orated honeycomb lattice.^ 
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FIG. 3. Phase diagrams for the decorated honeycomb lattice 
with t and t' real in the absence of a staggered on-site poten- 
tial and no Rashba coupling (from Ref. 84 ). Several filling 
fractions / are shown (lower left corner). For fixed / and 
Aso it is possible to drive a transition between a topological 
insulator and a non-topological phase by varying the ratio 
t'/t. 

We now show that the TI phase at 1/2 filling is 
"connected" to the chiral spin liquid phase of Q via 
the Kitaev-Haldane-Kane-Mele correspondence demon- 
strated on the honeycomb lattice by Lee, Zhang, and 
XiangP^ We do this by showing that the Hamiltonian 
describing the band structure of (|9| (with "normal" spin- 
less fermions replacing the Majorana fermions) with only 
purely imaginary first-neighbor hopping can be adiabat- 
ically deformed to one spin component of ( 10 1 (with 
"Hcdw = %R = 0) which has purely real first-neighbor 
hopping and purely imaginary second-neighbor hopping. 
This is done by applying a combination of gauge trans- 
formations and adiabatic tuning of the real and imagi- 
nary components of first and second neighbor hopping 
parameters. The procedure is shown in Fig. [4j and the 
panels in Fig. [5] explicitly show that the Chern number 
remains unchanged and the gap remains open in this pro- 
cess, establishing the topological equivalence, and the 
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FIG. 4. (color online) (From Ref.[84 .) Schematic illustra- 
tion of the continuous path which adiabatically connects the 
representative free fermion model of the Kitaev model on the 
decorated Honeycomb lattice with the spinless model with 
real t, t' and Aso- The patterns (o) and (i) are equivalent 
and correspond to Kitaev's model whereas (v) represents the 
model with real t, t' and Aso- The adiabatic deformation does 
not lead to a gap closing and the Chern number stays con- 
stant. See Fig. [5] This establishes the topological connection 
between the two models. 
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FIG. 5. (color online) Contour plot of (a) the Chern number 
and (b) the gap of the model as a function of {<j>, Aso) at half- 
filling for J = J'. Here <f> is the phase angle on the nearest 
neighbor hopping parameters given in Ref. [SI]. Also shown is 
a possible path which adiabatically connects the flux patterns 
(iii), (iv), and (v) defined in Fig|4] 

It should be clear from the examples of the decorated 
honeycomb lattice and the honeycomb lattice above that 
the Kitaev-Haldane-Kane-Mele correspondence is rather 
general. However, it is also true that one can invert 



the argument to use the existence of a TI in an s-band 
model on a particular lattice to infer a chiral spin liq- 
uid with a finite Chern number in a Kitaev-type model. 
Indeed, we have used the work of Guo and Franz that 
established a TI on the kagome lattice^ to motivate the 
study of a Kitaev-type model on the same lattice.^ Not 
only did we find a new exactly solvable chiral spin liquid 
model with finite Chern number, but also a phase that 
possess a stable spin Fermi surface.^ We note that un- 
stable spin Ferm i surfaces have been reported in other 
modelsPEMISa 

Because the kagome lattice has sites with four links 
(as opposed to three on the honeycomb, decorated hon- 
eycomb, and decorated square lattice), the effective spin 
in a Kitaev-type model must be at least 3/2p3 These 
higher-spin Kitaev models require that the spin-1/2 Pauli 
spin matricies be generalized to Gamma matriciesP*' The 
spin Hamiltonian we studjff^ on the kagome lattice is 

+j> A J2 r! 5 rf + J2 r ' 5r f' (ii) 

(ij)eA (ij)ev 



where we have distinguished the nearest neighbor cou- 
plings Jij as J a, Jv and JL 



as J ^ 



and based on 



whether the link (ij) belongs to an up (A) triangle or 
down (V) triangle (the kagome lattice can be viewed 
as a corner sharing network of up and down triangles), 
and (ij) is taken in the counterclockwise sense for each 
triangle. Locally, the five T-matricies satisfy a Clif- 
ford algebra, {Tf,T b t } = 2S ab , where a, b = 1,...,5, and 
r a& _ [r a ,r b ]/(2«). In terms of the components of the 
spin 5 = 3/2 operatorsj^EnH 

r 1 = -^{^,s*},r 2 = -^{s z ,s*},t 3 = ^{s x ,sy}, 
r 4 = ^[(s x f - r 5 = (s z ) 2 - j. (12) 



With the identification (12), it is clear the model (11) 



has a global Ising spin symmetry under 180° rotations 
about the z-axis, and possesses time-reversal symmetry 
(TRS) (although TRS will be spontaneously broken in 
the ground state from the triangular plaquettes as we 
described earlier) in addition to the translational and 
threefold rotational lattice symmetry mentioned above. 



In order to solve the higher-spin T-matrix model (11 ), 
a Majorana representation similar to that for spin-1/2, 
Q, is employed 



are similar to the spin-1/2 version and one obtains a Ma- 
jorana tight-binding model on the kagome with imagi- 
nary hopping parameters!-^ Solving the model on a strip 
geometry in the gapped phase with a finite Chern num- 
ber (±2) yields the results shown in Fig. [6 The gauge 
dependent but stable Fermi surface in a gap ess phase is 
shown in Fig. [7] 

(b) 
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FIG. 6. (a) Band structure on a cylindrical geometry for Ja = 
Ja = 1.0, Jv = Jy = 0.8, J5 = 0. There are two gapless 
chiral Majorana edge states (dotted lines) which overlap on 
each other, (b) For J A = 1.0, J' A = 0.6, J v = 0.9, Jv = 0.5, 
and J5 = 0.1, the two gapless edges states separate. These 
ground states are thus CSLs with a spectrum Chern number 
(±2) and the vortices obey Abelian statistics. From Ref. |18|. 





K 



(13) 



FIG. 7. The Fermi surface (solid line) for the 
flux configuration {|,|,7r} and { Ja, Jv, J'a, Jv> Js} = 
{1.0,0.3,0.8,0.5,1.4}. (See Ref.pl] for a description of no- 
tation.) The dashed hexagon is the Brillouin zone boundary. 
Note that there is only one Fermi pocket for this set of param- 
eters and the three pockets shown are related by a reciprocal 
lattice vector and are thus equivalent. 



with a — 1,2,3,4. There are thus 6 Majorana species 
on each site i: {£?• , , £f , £ 4 , Cj, d{\. The Majorana rep- 
resentation enlarges the spin-3/2 Hilbert space, so that 
one must enforce the constraint Di = —T\Y^TjTfT\ = 
—i£i£i£f£iCidi = 1, namely for any physical state |\t , ) p hys 
Di\ i S) p hys = |^)phys for any i. The remainder of the steps 



The results of this subsection have demonstrated 
the Kitaev-Haldane-Kane-Mele correspondence and illus- 
trated how TIs can be generated from a certain class 
of QSLs, and vice-a-versa. Thus, we have established 
a method of generating certain types of desired states. 
Moreover, along the way we have also discovered new 
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QSLs, such as the gapless varieties with a stable spin 
fermi surfaceP^In all these "non-interacting" models, the 
topological relationship between the TI and QSL models 
is based on the Chern number and an underlying band 
structure of free fermions. In the next subsection we turn 
to the study of interacting models in two dimensions. 

As a final remark, we note there are also deep 
connections between Kitaev models, the quantum 
Hall effect, top ological insula tors, and topological 
superconductivity.EMllIIinSHinZl The topological Kitaev 
models 19 and the TI insulators (including the integer 
quantum Hall effect ^pHHUl s how that disorder can play a 
useful role in stabilizing topologically ordered states and 
could play a role in guiding the experimental discovery of 
topological quantum spin liquids !^ The non-interacting 
models discussed in this section can be viewed as mean- 
field states in a more general interacting model. 



C. Interacting lattice models 

While the Kitaev models are interacting in terms of 
the underlying electronic degrees of freedom (local mo- 
ments), we discussed them under the heading of "non- 
interacting" lattice models because they admit a free 
(Majorana) fermion representation and therefore can be 
understood in terms of the properties of their band struc- 
ture. There are few exactly solvable interacting models 
in two dimensio ns, a nd those that exist often possess spe- 
cial symmetries?^ Thus, one must look for approximate 
methods to address the physics of interactions. 

For a certain class of Hamiltonians (those without a 
fermion sign problem), quantum Monte Carlo is a pow- 
erful theoretical method. Given the central role the hon- 
eycomb lattice has played in the understanding of TIs 
and QSLs, it is natural to begin an investigation of inter- 
acting models on this lattice. In the last year the simplest 
interacting model, the Hubbard model (no spin-orbit 
coupling terms), has been shown to possess a (appar- 
ently gapped) sp in-li quid phase at half-filling and inter- 
mediate couplingP^H Moreover, this spin liquid has been 
shown to be robust in the presence of a seco nd-n eighbor 
spin-orbit couplin g in quantum Monte CarlcW^* and dy- 
namical mean- fielcP^ studies of the so-called Kane-Mele- 
Hubbard model. Related quantum Monte Carlo studies 
on bulk and edge instabili ties in t he Kane-Mele-Hubbard 
model have been reportedp^^as well as closely related 
exact diagonalization studies in spinless systems wit h 
nearest-neighbor interactions for small system sizes 

Aside from bosonization studies of interactions in 
the one-dimensional helical liquids that ex ist along the 
boundary of a two-dimensional Tlpfiffilal most studies 
of interacti ons in bulk hav e relied on some type of mean- 
field theorypmMEMIHI] 0ne f the most important re- 
sults to emerge from these studies is that topological 
insulators can be spontaneously generated from interac- 
tions at the mean-field level even when there is no micro- 
scopic spin-orbit coupling present. This was first explic- 



itly shown by Raghu et al. on the honeycomb lattice.^ 
In model without any intrinsic spin-orbit coupling, TI 
phases are produced by the spontaneous generation of 
spin-orbit coupling from first and seco nd-neighbor in ter- 
actions treated at the mean-field levelP* 7 ^ * 86 * 1241 ^^ 

The stability of a non-interacting band structure to 
interactions depends to a large extent on the form of 
the low-energy dispersions and band crossing points. For 
example, it is well-known that at 1/2 filling on the 
honeycom b la ttice the non-interacting dispersions are 
Dirac-likeP ^* Dirac po ints are perturbatively stable to 
interactions! 17 * 86 * 88 * 124 ^ In practice, this implies that a 
critical interaction strength is required to drive a phase 
transition away from the gapless stated This is quali- 
tatively consistent with both mean -fiel d studies^ 7 * 124 * 129 * 
and quantum Monte Carlo studies .1221 By contrast, two 
bands that touch quadratically or flatter are in finites- 
imally unstable to interactions in two dimensions ! 86 * 12 ^ 
Quadratic band touching points are known to occur 
on many of the lattices that support TIs in sim- 
ple s-band models: decorated honeycomb,^ kagomep^ 
checkerboard,^ square-octagonj^] and rubyP^ As we 
show below, the cases of quadratic band touchings points 
often have topological phases a s th e leading instabili- 
ties with respect to interactions This effectively en- 
larges the candidate systems for TI physics to include 
those without strong spin-orbit coupling but with "flat" 
band touching points in the non-interacting band struc- 
ture. An othe r notable example is bilayer and multilayer 
grapheneP^2] 

The kagome lattice serves as an important exam- 
ple case. Its non-interacting band structure in an s- 
band nearest neighbor hopping model possess both Dirac 
points (1/3 filling) and quadratic band touching poin ts 
(2/3 filling) P We study the following HamiltoniarP^ 

"Hkagomc = —t c\ a Cj a + U ^ U^n^ 
(id) 

+Vl ^ 71 i n 3 + riiUj + V 3 ^ n * n 3 ( 14 ) 

(id) ((id)) (((id))) 

at the mean-field level. Here, c\J annihilates (creates) 
a fermion on site i with spin <r =j, J,, n: la = c\ a c ia and 
rii = ^2i a n ia- The sums run over nearest-neighbor 
second-neighbor or third-neighbor bonds (((i,j)))- 

The hopping amplitude is denoted by t and the parame- 
ters Vi , V2 , and V3 quantify the nearest-neighbor, second- 
neighbor and third-neighbor repulsion, respectively. 

To perform the mean-field ca lcul ation, we decouple the 
on-site interaction according td^* 

riifnn « n lt (n 4 ) + (ntf)mi - (n lT )(n. 4 ) 
" c !t c 4( c 4 c iT> " (4t c '4) c 4 c *T + (4r c 4)( c l c *T>-( 15 ) 

We assume the mean-field solutions are described by a 
colinear spin alignment and therefore, without loss of 
generality, we set (cLc 4 ) = (ctc^) = 0.^231 The further- 
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neighbor interaction is decoupled in a similar wayP^iJ 

mrij w m{rij) + (ni)nj- - (n 4 )(n 3 ) - ^ (cj Q c^{c] p c ia ) 

a/3 

+ (4 a c jP) c ]p c ic* ~ ( c L c j/3>( c ]/3 C »)) • ( 16 ) 

The results of the mean-field calculations are shown in 
Fig. [8] and [9] For a detailed description of the phases 
appearing in the phase diagrams, see Ref. |124j . 
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FIG. 8. Filling fraction 1/3 on the kagome lattice with Dirac 
points in the non-interacting band structure. The U-V2 phase 
diagram for Vi = and V3 = 0.4i in \14\ . An interaction- 
driven TI appears for finite U and V2 and requires a critical 
interaction strength to appear, and some fine-tuning of pa- 
rameters. Solid lines indicate first order and dashed lines 
second order transitions. From Ref. 1 2 I . 
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FIG. 9. The U - Vi phase diagram for V 2 = V 3 = at 2/3 
filling fraction (quadratic band touching point) in \14) . Note 
that for small generic interactions the leading instability is a 
topological phase, either TI or the quantum anomalous Hall 
phase (QAH). From Ref. [124]. 

The important message we would like to emphasize is 
that for Dirac points, Fig [8j topological phases require a 
critical interaction strength that can be rather la rge, and 
appear to require some amount of "fine tuning" P^ For 
example, note that V\ = in Fig. [8] On the other hand, 
when quadratic band touching points are present in the 



non- interacting band structure, Fig. [9j topological phases 
can appear as the leading instabilities, without any "fine 
tuning" of the interactions. Thus, two-dimensional band 
structures with flat band touching points but small (or 
even vanishing) spin-orbit c oup ling are excellent candi- 
dates for topological phasesP^ General arguments and 
functional renormalization group studies (that treat fluc- 
tuations) suggest the mean-field results may be reliable 
and that the TI phase is likely preferred over the quantum 
anomalous Hall (QAH) state when they are degenerate 
at the mean-field levelP^ 

A final point on which we would like to remark con- 
cerning interacting topological insulators in two dimen- 
sions is how interactions may help identify the phase it- 
self. For example, if two edges of different quantum spin 
Hall systems are brought close to one another the edges 
will interact via the Coulomb interactions of the elec- 
trons in each edge. Under some circumsta nces , this can 
lead to novel quasi-one dimensional phasesP^H However, 




I, 



FIG. 10. (color online) Schematic of a drag measurement 
between two QSH systems from Ref. |12l) . A current Ji is 
driven along the upper edge of the lower QSH system and 
through electron-electron interactions a voltage V2 is induced 
in the lower edge of the upper QSH system. A magnetic field 
h is applied in the plane of wires, perpendicular to the spin 
quantization axis (assumed perpendicular to the plane of QSH 
systems). Time-reversed Kramer's pairs are indicated for the 
two edges. A QSH on top of QSH geometry could also be 
used. 



if the interactions are not too strong (the edges are far 
enough apart), a current driven in one edge will result 
in a voltage drop alo ng the o ther edge in the so-called 
"Coulomb drag" effect The properties of the edge 

are reflected in the temperature and magnetic field de- 
pendence of t he dra g. For identical edges described by 
a helical liquid, 68 69 the Coulomb drag has a novel mag- 
netic field dependence that is "turned on" in the pres ence 
of a time-reversal symmetry breaking magnetic field P^J 



rrj oc h T 



4K-3 



(17) 



where h is the strength of the applied magnetic field 
and K is the Luttinger parameter of the helical edge 
systemP^H The role of the finite magnetic field is to open 
a back-scattering channel for electrons. There is als o a 
weak field dependence of the Luttinger parameter K)r^ 
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D. 



Fractional topological insulators and flat band 
fractional quantum Hall effect 



The discussion in Sec. |III B| which related the integer 
quantum Hall effect to TIs immediately raises the ques- 
tion of whether there is a "fractional" analog to TIs that 
depends in an intrinsic way on electron-electron interac- 
tions in the same way the fractional quantum Hall ef- 
fect does. In particular, is there a time-reversal invariant 
topological insulator that is not adiabatically connected 
to a free electron model? Theoretically, the answer is 
affirmative!^^ 

In two dimensions, we are guided in a study of frac- 
tional topological insulators by the fractional quantum 
Hall effect. If one has in hand a spinless fermion model 
that realizes a fractional quantum Hall effect, an S z con- 
serving, time-reversal invariant model can be straightfor- 
wardly constructed by generalizing the Haldane-Kanc- 
Mele correspondence to the interacting caseP2B3l Only 
very recently has the explicit construction been under- 
stood for the non-S" 2 conserving case (in the presence of 
Rashba spin-orbit coupling, for example).^ 

A key to understanding lattice models of fractional 
topological insulators in two dimensions is thus un- 
derstanding the fractional quantum Hall effect in lat- 
tice models. (While this appears necessary, it is 
not sufficien t for con structing a general many-body 
wavefunctionP^l^ShliUj p or a w hile, this understanding 
eluded the community whose picture was based on trial 
wavefunctions and effective low-energy theories,^ but in 
the past year important progress has been made. 

Attention has focused on a class of insulators with 
nearly flat bands that possess finite Chcrn numbers. 
With a finite Chern number, a partial filling of the flat 
bands can lead to a fractional quantum Hall effect be- 
cause inte raction energy will always dominate kinetic 
energyfSH To date, only a few lattice models have been 
proposed which are expec ted to lead to a fractional quan- 
tum Hall effect \MlMMil The relevant figure of merit in 
such models is the ratio of the band gap to the bandwidth 
of the flat band with a finite Chern number. In the model 
we discuss below, we find this ratio can be as high as 70, 
which is among th e largest i n the models reported in the 
literature thus farPESMHI] 

We study^a simple s-band tight-binding model on the 
ruby lattice shown in Fig. |11| For fermions with spin, the 
Hamiltonian is given by H = H ru by + ^ruby-so, where 



H 



ruby 



-t 



E 

<JeA,cr 



C ia C 3<? 



E 

A->-A,, 



and 



■H ru by-SO = ^2 E "<J ' S al3 C L c ] + »*3 E ^J S a/3 c L c j/3 

(19) 

on the ruby lattice shown in Fig. [ill Here <y ia jci a is 
the creation/annihilation operator of an electron on site 
% with spin a. As indicated in Fig. |ll[a,b), t and t\ 
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FIG. 11. Schematic of the ruby lattice and illustration of 
the nearest-neighbor hopping, t,t\ (real) and (complex), 
and the three types of "second-neighbor" spin-orbit coupling 
or hopping indicated by the dashed or dot dashed lines, t2,t^ 
and t± r (from Ref. |85|). (a) The spin-orbit coupling strength 
ti within a hexagon, (b) The spin-orbit coupling strength t$ 
within a pentagon composed of one triangle and one square. 
The hoppings t^ and t$ axe present on all such bonds of the 
type shown that are consistent with the symmetry of the lat- 
tice, (c) Schematic of the hopping parameters used to obtain 
a flat band with a finite Chern number and W/E g ~ 70. 



are real first-neighbor hopping parameters, and t2, t% 
are real second-neighbor hoppings (these appear with the 
imaginary number i in Eq.(19) making the total second- 



neighbor hopping purely imaginary and time-reversal 
symmetric). The quantity i/y is equal to 1 if the elec- 
tron makes a left turn on the lattice links during the 
second- neighbor hopping, and is equal to -1 if the electron 
make a right turn during that process. As is clear from 
FigjTTj the unit cell of the lattice contains six sites, so six 
two- fold degenerate bands will result. In addition to the 
real hopping parameters i, ti in (18 1, symmetry also al- 



lows complex, spin-dependent nearest neighbor hopping 
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with imaginary components t'jii, as shown in Fig 
These parameters lead to a variety of time-reversal in 
variant phases (including TI) and a rather complex phase 
diagram.^ 

Here we are mainly interested in the spinless case with 
broken time-reversal symmetry and the hopping param- 
eters described in Fig. 1 1 c) . For spinless fermions we 



(18) study (18) with 



t' = t + ia z ti, 
t'i = t\ r + ia z tu, 



(20) 



and add the hopping inside the square in the diagonal 
directions, labelled t± r . This is shown schematically in 
Fig. |ll[ c) . A nonzero Chern number and nearly flat band 
occurs in this case. For example, with ti — 1.2t, t\ r = 
—1.2t,t\i = 2.6t,t 4r — —1.2t the Chern number is -1 and 
the gap = 2.398i, and band width = 0.037£, which gives 
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FIG. 12. (a) The energy bands without any spin-orbit cou- 
pling and with it tuned on, given by Eq.( |18[ l with t = t' 
real along high symmetry directions. Note that there is a 
Dirac point at K for 1/6 and 2/3 filling, and a quadratic 
band touching point at V for 1/2 and 5/6 filling. (The 
underlying lattice is triangular, as it is for the honeycomb 
lattice.) Note also the flat bands along the Y — M direc- 
tion at 1/2 and 5/6 filling, (b) The energy bands with 
U = 1.2t,t lr = -1.2t,tu = 2.6t,t 4 r = —1.2* in Q. The 
lowest energy band has Chern number -1. The ratio of the 
band- width to band gap is approximately 70 (from Ref . \85\ ) . 



W/Eg « 70, and therefore makes an excellent candidate 
to realize a fractional topological insulatorP3 The corre- 
sponding band structure is shown in Fig. 12 b). The spin- 
less band structure for (18) with t' = t = 1 is shown in 



Fig. [I2[a). Note the interesting band crossing/touching 
points and flat portions of the bands even in this case. 

A new direction has emerged in the search for sys- 
tems with flat bands with conditions favorable to re- 
alizing fractional topological insulators or fractional 
quantum anomalous Hall states-transition metal oxide 
interfaces.^ Research in this direction is rapidly evolv- 
ing with a variety of interesting interaction-driven phases 
possiblePHU 

As rich as the physics of topological insulators and 
quantum spin liquids is in two dimensions, the behav- 
ior is even richer in three dimensions. This is be- 
cause the three-dimensional time-reversal invariant topo- 
logical insulators are qualitatively different from their 
two-dimensional counter parts, the latter of which can 
always be viewed as being derived from an underly- 
ing quantum Hall effect. The so-called "strong topo- 
logical" insulator d 1 ' (non-interacting) are an entirely 



distinct t hree-dimens ional phase with novel response 
properties! 8 ' 51 ' 143 ^ 4 ^ Moreover, once interactions are 
brought into the picture the physics is richer still. We 
now turn to three dimensions. 



IV. THREE-DIMENSIONAL SYSTEMS 

One of the salient features of TIs is that their bound- 
aries possess a topologicall y protected "meta llic" state 
that is robust to disorderPH' 69 ' 143 ' 149 ' ^ Both the 
two dimensionaf^H^and the three dimensiona l 8 ' 1441 * 1 ^ 
boundaries have been shown to exhibit interesting re- 
sponses to perturbations. In the experimental litera- 
ture on TIs, the three dimensional systems have domi- 
nated. The weakly-interacting nature of topological in- 
sulators has enabled accurate pre dictions based on den- 
sity functional theory^l 147 * 152 * 15 -^for a wi de range of two 
and three dimensional systems[ 138 ' 154 ' E^and experiment 
has followed with confirming data in a lar ge and rapidly 
growing number of instances.' 79 ' 80 ' 168 -'* 1 ^ 7 -' There are also 
theoretical predictions for systems closely related to TIs 
with Weyl fer mions in their bulk and unusual bound- 
ary excita tion s,' 32 ' 178 ' * 1 ^' as well as crystalline topological 
insulators! 1 ^ 3 -' 

In this section we will describe systems that are too 
strongly interacting to allow a naive application of den- 
sity functional theory. On the most exotic side, three- 
dimensional versions of exotic fractional topological in- 
sulators with a non-tr ivial ground state degeneracy have 
been proposedP^E^ For interacting systems treated at 
the mean-field level, a spontaneous generation of spin- 
orbit coupling (with topologic al de fects possible) can oc- 
cur, similar to two dimensions 

Here, we will focus on semi-realistic models 
that are motivated by the physics of transition 
metal oxides, famous in re cent y ears for high- 
temperature su perco nductivit y 188 1 189 1 and colossal 
magnetoresistanceP^ni These materials typically involve 
3d orbitals on the transition metals and lead to strong 
Coulomb interactions due to their localized nature 
(compared to 4c? and 5c? orbitals). However, Ad and 
5c? orbitals in layered perovskites such as Sr 2 Ru04, 
Sr 2 Rh0 4 , Sr 2 Ir0 4 , Na 2 Ir0 3 , and the hyperkagome 
Na^^Os are more spatially extended and thus the 
Coulomb interac tion is typically weaker than those 
with 3d orbitalsP^H The hyperkagome Na^^Os is 
an importan t th ree-dimensional quantum spin liquid 
candidate™!! 

The more extended nature of the 4c? and 5c? orbitals 
compared to the 3c? orbitals leads to a greater level split- 
ting in a crystal field, and enhances their sensitivity to 
lattice distortions! 1 ^ 7 -' In many oxides, the transition ions 
are surrounded by an octahedron of oxygen atoms, M06, 
where M represents a transition metal ion. The crystal 
field splits the 5 degenerate (neglecting spin for the mo- 
ment) G?-orbitals into two manifolds (see FigJT3fc) : a lower 



lying t 2& 



(d X y,dy Z ,d zx ) manifold and a higher lying 
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FIG. 13. (Color online) From Ref.^Tj. (a) An illustration 
of the pyrochlore lattice which is composed of corner sharing 
tetrahedra. Transition elements are indicated by black solid 
circles, (b) Each transition ion is surrounded by an oxygen 
octahedron shown by six solid blue (dark grey) circles. A 
transition ion is located at the origin of the local coordinate 
and is shown in black. We study a trigonal distortion pre- 
serving C3 symmetry applied along the [111] direction (or its 
equivalent), shown by two yellow (grey) faces, and an elon- 
gation preserving C4 symmetry along the z-axis of the local 
coordinate, (c) A schematic representation of the splitting of 
the bare atomic d-levels (1), due to a cubic crystal field arising 
from the octahedral environment (2), unquenched spin-orbit 
coupling in the t2 g manifold (3), and a distortion of the oc- 
tahedron (4). The values of the splittings in (4) depend on A 
and A3 4. 



(d,3 Z 2_ r 2,d x 2_ y 2) manifold! 198 * 1 " ! The energy separation 
between the t^g and e g levels is conventionally denoted 
"lODq" and is typically on the order of ~l-4 eV, which 
is large compared to many 3d compounds I 2221 

Besides the crystal field, the relativistic spin-orbit cou- 
pling is another energy scale that results from the large 
atomic numbers of heavy transition elements. While in 
the absence of spin-orbit coupling the on-site Coulomb in- 
teraction is of the same order as the band width 12211202] j n _ 
elusion of str ong spin-orbit coupling modifies the relative 
energy scalesPSl Thus, for materials with 4c? and partic- 
ularly 5d electrons, one expects the appearance of novel 
phases with unconventional electronic structure and or- 
der due to the characteristic energy of spin-orb it couplin g 
approaching that of the Coulomb interactions! 27 * 203 * 204 * 

In this section we focus on the interplay and com- 
petition between strong correlation effects, spin-orbit 
coupling, and lattice distortion that is expected to be 
important in heavy transition metal oxides with a py- 
rochlore lattice! 2 ^* In the heavy transiti on meta l oxides 
one expects both the spi n-orbit coupling 205 2 ® and the 
lattice distortion energies^ 197 * 207 ! to be of the order of 
0.05 — 0.5 eV, while the interaction energy is typically at 



the higher en d of this scale to somewhat larger, 0.5 — 2 
eV! 178 * 179 * 205 * While the phase diagram of an interacting 
undistorted pyro chlore model with j = 1/2 has already 
been studiedpSSl W e expand those results to include the 
effects of distortions of the local octahedra on the phase 
diagramP3 

We also investigate pyrochlore oxides at different d- 
lcvcl fillings with the Fermi energy lying in the quadru- 
plet j — 3/2 manifold, which has not been considered in 
previous works One of our motivations is to see if the 
j = 3/2 manifold can also re alize the interesting Mott 
phases of the j = 1/2 manifold P^* We find that, indeed, 
these exotic phases can be realized for the j = 3/2 man- 
ifold. Moreover, we find that for the j = 1/2 manifold 
"weak" topological variants of the exotic Mott phases can 
also appear in the phase diagram when certain types of 
lattice distortion are present. 

To study the effects of lattice deformations 112Z12Q8J we 
assume that the octahedron surrounding an ion can be 
distorted in two ways: (1) a trigonal distortion preserving 
local C3 symmetry and (2) an elongation (expansion) of 
octahedra preserving local C4 symmetry. (See Figp^jp.) 
The former has been argued to be rather common and 
can be described by th e fol lowing Hamiltonian on each 
transition metal ion siteP^ 



tit 



A 3 (d\ z d zx + d\ z d xy + d\ x d xy ) + h.c, (21) 



where A3 parameterizes the strength and sign of the C3 
preserving distortion, a nd t he C4 elongation/contraction 
splitting is described by^' 



Hci = A 4 ^ 2 = A 4 (n yz + n zx ), 



(22) 



where A4 parametrizes the strength and sign of the dis- 
tortion, and l z is the z component of the effective angular 
momentum of the ti g orbitals related to the occupation 
of the d xy orbital by n xy = n c i — (l z ) 2 w hich follows from 



the constraint rid — n xy 



EHl For an elongation 
of the tetrahedron, A 4 < 0, and for a compression of the 
tetrahedron, A 4 > 0. Trigonal distortions appear to be 
more common in real materials, and the magnitude of 
the energy splittings can be crudely estimated from den- 
sity functional theory calculations based on X-ray deter- 
mined positions of oxygen atoms around the transition 
metals. We are not aware of detailed calculations of this 
type for the id and 5g? pyrochlore oxides, but closely re- 
lated 3d syste ms ap pear to have splittings on the level 
of 0.01 — 0.5 eVpSZl We take this as crude estimate, with 
the larger end of the energy scale probably more likely 
for the more extended 4d and 5d orbitals. 

Thus, the local Hamiltonian describing the t 2g orbitals 
on each site is 



u 



local 



Ti-so + T~Ltr. 



H 



el 3 



with 



T~Lso 



-XI ■ s, 



(23) 



(24) 
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where I = 1 and s = 1/2 describe the orbital and spin 
degrees of freedom, and A > parameterizes the strength 
of the spin-orbit coupling. That the t 2g orbitals can be 
effectively described by angular momentum / = 1 comes 
from the projection of the rf-orbit al angu lar momentum 
into the local basis of t 2g manifoldP^ESS] 

The Hamiltonian (23) can be easily diagonalized and 



its eigenvectors describe a projection onto the spin-orbit 
plus distortion basis. We will denote the projection by 
a matrix M, which contains all the information about 
the spin-orbit coupling and the distortion of the octa- 
hedra (all assumed identical so translational invariance 
is preserved). Moreover, due to the presence of time- 
reversal symmetry, the eigenvectors form Kramers pairs. 
A schematic representation of splitting t 2Q upon includ- 



ing the terms in Eq.(23) is shown in Fig 13:. 

We now turn to a derivation of the effective Hamilto- 
nian. We first assume A = A3 = A4 = 0, i.e. neglect 



the contributions in ( 23 1 . To obtain the kinetic terms 



of the Hamiltonian, we need to describe the t 2g orbitals 
of a single ion in the local coordinate system defined by 
the octahedron of oxygen atoms surrounded the ion, and 
we nee d th e p-orbitals of oxygen in the global coordinate 
systemPS^l The hopping of electrons from one transition 
metal ion to a nearest-neighbor transition metal ion is 
mediated by the oxygen p-orbitals. (We note that for the 
relatively extended 5d orbitals direct overlap ma y als o be 
important, as well as further neighbor hoppingPSI) We 
thus compute the p-d overlaps to determine the hopping 
matrix elements. The local and glo bal axes are related 
by a set of rotation matrices! 197 * 203 ^ The combination of 
rotation matrices and d-p overlaps gives rise to the fol- 
lowing Hamiltonian: 



T~Ld — £d 



t 

i'ya' 



<27fT,i'7 / er' > 



T U ' , ,d\ di / / 



(25) 

where i, 7, and a in the sums run over lattice sites, t 2g 
orbitals (xy, yz, zx), and spin degrees of freedom, respec- 
tively. The £d stands for the onsite energy of the degen- 

v 2 

erate t 2ci orbitals, and t — v ^ is the unrotated hopping 
amplitude depending on the overlap integral V v d-n and the 
energy difference between p and d orbitals. The parame- 
ter t sets the basic hopping energy scale in the problem. 
Without loss of generality we set Ed = 0. 

The effect of spin-orbit coupling and distortion are in- 
cluded via the projection of Hamiltonian in Eq.(25) into 
the eigenvectors of the local Hamiltonian in Eq. ( 23 ) using 
matrix M as follows: 



ia 



<2Q,i' a' > 



/C ia c i'at' 



(26) 



where v a (a — 1, 6) stands for the six eigenvalues 
of local Hamiltonian (23), and the matrix T describes 



the hopping between sites given in the local basis via 
r = M*TM T . The c\ a {c ia ) is the creation (annihilation) 
operator of an electron at site i and in local state a. 
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FIG. 14. Phase diagram of the j = 1/2-band model cor- 
responding to rid = 5 with positive A3 = 2t (upper panel) 
and negative trigonal distortion A3 = — 2t (lower panel). The 
dashed line separates the rotor condensed phases (below) from 
the uncondensed phases (above). We set t = 1, and the phases 
labeled are as follows: Strong topological insulator (STI), 
Weak topological insulator (WTI), Gappless Mott insulator 
(GMI), Topological Mott insulator (TMI), Weak topological 
Mott insulator (WTMI) and Metallic phases. From Ref.[27]. 



Finally, we add a Coulomb interaction to obtain 

U 
2" 



where U is the on-site Coulomb interaction and rid is 
the number of electrons on the 5d orbital of the transi- 
tion metal ion. In the remainder of this section, we in- 
vestigate the zero-temperature phase diagram of the full 
Hamiltonian (27), which includes the spin-orbit coupling 
and lattice distortions in (23). A full description of our 



results is given in Ref . [2"T] . 

We apply the slave- rotor mean-field theory developed 
by Florens and Georges^^fil ^ treat the effect of weak to 
intermediate strength Coulomb interactions in the regime 
where the charge fluctuations remain important. In this 
theory each electron operator is represented in terms of a 
collective phase, conjugate to charge, called a rotor and 
an auxiliary fermion called a spinon as 



H>i 



fit 



(28) 



where Ci a is the electron destruction operator at site 



with quantum number a, representing the states in (27) 
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FIG. 15. Phase diagram of the j — 3/2 band model with 
rid = 3, including the trigonal distortion of the octahedra. 
The labeling of the phases is the same as that used in Fig |14| 
In the upper panel A3 = 2t. In the lower panel the interac- 
tion is fixed at U — At and the strength A3 of the trigonal 
distortion is varied, illustrating possible phases that may arise 
upon the application of pressure to a real system. All ener- 
gies are expressed in units of t. In both phase diagrams the 
dashed line separates the rotor uncondensed phase (above) 
from the condensed phase (below). We note the "pocket" of 
STI around A w 1, A3 « 1.5 in the lower figure has a numeri- 
cally difficult to determine boundary with the metallic phase; 
we have present our best assessment. From Ref.|27|. 



The factor e lBi acting on the charge sector is a rotor low- 
ering operator (with 0j a bosonic field), and fi a is the 
fermionic spinon operator. The product of the two results 
in an object with fermi statistics, needed for the electron. 
Note the rotor part only carries the charge degree of free- 
dom while the spinon part carries the remaining degrees 
of freedom a. Therefore, an electron has natural spin- 
charge separation if a is spin in this representation. Sub- 
stituting (28 1 into (27 1 and applying mean-field t heory 

andli~5l 2 ' 7 
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results in the phase diagrams shown in Figs 
The two to three dimensional crossover and issues asso- 
ciated with gauge fluctuations beyond mean-field theory 
have been discussed in Ref. 210 . The method correctly 
recovers the U — > limit in mean-field theory. 

There are several features that are important to em- 
phasize. First, we have identified a new phase-the weak 
topological Mott insulator (WTMI)-in the regions shown 
in Fig. 14'b)p3This phase is expected to have gapless 
thermal transport without charge transport along a cer- 
tain class of topological defects, as described in Ref. [35J . 



Also, we expect a rather rich phenomenology in the spin 
sector of WTMIs, similar to that recently emphasized in 
weak topological insulators (WTI)PSl Second, it is clear 
from Fig.[l4]that the sign of the trigonal distortion makes 
a big difference in the complexity of the phase diagram, 
with the negative trigonal distortion producing a more 
interesting phase diagram. Third, comparing Figs. |14| 
and 15 one sees that the j = 1/2 and j = 3/2 mani- 
folds behave differently for the same distortion. Fourth, 
the lower panel of Fig. [15] shows the phase diagram in 
the distortion-spin-orbit coupling plane. The distortion 
axis can be viewed as a pressure axis. Thus, the applica- 
tion of pressure to transitional metal oxides can drive the 
system through a complex set of phases, some of which 
are topological. This is potentially a useful experimental 
parameter when searching for some of the exotic phases 
predicted here. 

Our slave rotor results make a natural connection to 
quantum spin liquids in three dimensions described by 
fermionic excitations: All the phases above the dashed 
lines in Figs \14\ and \ lh\ are time-reversal invariant spin 
liquids. These phases have fully gapped charge degrees 
of freedom and spin degrees of freedom described by the 
fermionic spinons "/«*" in (28). At the mean-field level, 



the spinons have an associated band structure that can 
either be gapped and topological (like a TI) or gapless 
(like a metal). Thus, the spin sector has an intrinsic 
single-particle-like nature imposed on it by the formal- 
ism. Relaxing this assumption is one of the ingredients 
needed to obtain the most exotic fraction alized topologi- 
cal insulators in three dimensions JSIEHZI Indeed, a fruit- 
ful approach to studying and/or discovering exotic frac- 
tionalized topological insulators in three dimensions may 
be to deepen our understanding of quantum spin liquids 
in three dimensions, particularly those whose spinons are 
strongly correlate d and not well described within a single- 
particle formalism!^] In a gapless system, these would be 
"algebraic spin liquids" P One could then "add back" the 
charge degrees of freedom to obtain an exotic topolog- 
ical insulatorPS However, it should be emphasized that 
correspondences between TIs and QSLs described in this 
paper have all been between QSLs with excitations de- 
scribed by fermions. At the moment, it remains unclear 
how to relate spin liquids described by bosonic excita- 
tions (Schwinger bosons) to the topological insula tors , 
although some interesting parallels have been notedPi2l 



V. CONCLUSIONS 

In this work, we have emphasized connections between 
topological insulators and quantum spin liquids from a 
variety of viewpoints, with an emphasis on their com- 
mon topological structure. We have described results 
from one, two, and three dimensional systems. In one- 
dimension we have emphasized the important role that 
quantum entanglement, particularly the entanglement 
spectrum, can play in studying and classifying phases. 
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In two dimensions, we have demonstrated the Kitaev- 
Haldane-Kane-Mele correspondence and shown how it 
can play a useful role in finding topological insulators and 
interesting quantum spin liquids. In interacting systems 
we have described the different roles Dirac points and 
quadratic band crossing points play in the instability to 
spontaneously generated topological phases, even when 
there is no intrinsic spin-orbit coupling. We have given 
an example of a lattice model of spinless fermions that 
possesses a very flat band with a finite Chern number. 
This model and related ones can be used to help guide 
the discovery of fractional topological insulators in two 
dimensions. We have also emphasized transition metal 
oxide interfaces as a promising experimental system for 
the realization of interaction-driven topological phases. 

In three dimensions, we have focused on semi-realistic 
models appropriate for transition metal oxides with Ad 
and 5d electrons. We have studied the interplay of distor- 
tion, interactions, and spin-orbit coupling and presented 
rich phase diagrams for different d-shell fillings. We have 
also found a new phase, the "weak topological Mott insu- 
lator" , which is expected to reveal itself via topologically 
protected gapless modes along certain classes of defects.^ 

Looking ahead, it is clear that two of the chief 
challenges with three-dimensional TIs (which are most 
likely to appear in applications) are sample quality 
and reliable tunabilityP^ To date, the most pressing is- 
sue is to obtain samples that have a highly insulat- 
ing bulk (so that most electrical conductance is along 
the boundary). In recent months importan t adva nces 
have been made o n this fr ont. Both Bi2Te2S c! 213 * 214 l and 
Bii, 5 Sbo.5Tei.7Sei. j 215 | 216 | are reported to have more than 



70% of their conductivity on the surface, and there are 
indications that T I "w ork horses" like Bi2Se3 are improv- 
ing in quality, tooP^l 

In the regime of very strong electron interactions (be- 
yond where the slave-rotor formalism should apply), 
magnetic phases with broken time-reversal symmetry be- 
come more likelyP^ When spin-orbit coupling is strong 
in this regime, exotic, co nventionally o rdered phases 
are also expected to appea r j 195 l 196 l 204 l 218 l including spin- 
orbital liquid states E. 9 I 220 I It seems likely that the new- 
found capability to engineer gauge fields in cold-atom 
systems w ill play a role in our understanding of TIs and 
QSLs, too ' 221 " 226 ! Many novel behaviors certainly await 
identification, and keeping in mind the "connections" be- 
tween various topological phases seems a good strategy 
to help guide discovery. 
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